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Light vector-like quarks with non-renormalizable couplings to the Higgs are a common feature of 
models trying to address the electroweak (EW) hierarchy problem by treating the Higgs as a pseudo- 
goldstone boson of a global (approximate) symmetry. We systematically investigate the implications 
of the leading dimension five operators on Higgs phenomenology in presence of dynamical up- 
and down-type weak singlet as well as weak doublet vector-like quarks. After taking into account 
constraints from precision EW and flavour observables we show that contrary to the renormalizable 
models, significant modifications of Higgs properties are still possible and could shed light on the role 
of vector-like quarks in solutions to the EW hierarchy problem. We also briefly discuss implications 
fT^ of higher dimensional operators for direct vector-like quark searches at the LHC. 

o 

O , I. INTRODUCTION 

< 

II ever more stringent constraints from measurements of flavor, EW and Higgs observables are starting to directly 

CIh probe models addressing the naturalness problem of the SM. At present these experimental null-results are not yet 

^ conclusive, but viable new physics (NP) models' parameter space is becoming significantly reduced. In light of this, it 

has become crucial to focus on the minimal light particle content required to fulfil the naturalness conditions (see [2] 

for examples in supersymmetric theories). 

—^ Vector-like quarks are expected to be the lightest new degrees of freedom in models addressing the EW hierarchy 

^ problem by treating the light Higgs as a pseudo-goldstone boson of a global symmetry, broken explicitly by the SM 

gauging and the top Yukawa [3, 4]. In such models, the dominant quadratic divergences in one-loop corrections to 

' -j the Higgs boson mass coming from the top quark loop are canceled by contributions of dimension five operators of 

—^ the form H^ HQQ, where H is the Higgs doublet and Q a vector-like quark weak multiplet. In general, parametrizing 

• single- and double-Higgs interactions of an arbitrary number of quark flavors / in the mass eigenbasis 
■^ 

Cn rf =(^-y„■/^ + a;„■^)7l/^ + h.c., (1) 

.5^ where v = {\/2Gf)~^^^ — 246 GeV is the EW condensate, the condition for cancelation of one-loop quadratic 
S^ divergences from such interactions can be put into the following simple form 
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The recent LHC discovery of a Higgs boson [1] seems to experimentally complete the standard model (SM) picture 
of fundamental interactions. Although the SM works extremely well phenomenologically, the electro-weak (EW) 
hierarchy problem, as exemplifled by the extreme UV sensitivity of the Higgs potential, provides a strong theoretical 
motivation for contemplating beyond SM physics at the TeV scale. 

On the other hand, direct searches for on-shell production of new degrees of freedom at the LHC together with 



Y^^ix.,)-^ = j2\y^^\'- (2) 
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The appearance of new operator contributions already at dimension five is a singular feature of effective theories with 
new light vector-like fermions and intrinsically connected to the resolution of the SM hierarchy problem within such 
scenarios. These dimension five contributions are thus expected to represent the dominant probe of new dynamics in 
the UV. 
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Treating the Higgs as a composite field of a strongly interacting theory leads to the appearance of a number of 
dimension six operators affecting flavor, EW and Higgs observables (c.f. [5]). However, in models where the new vector- 
like quarks (possibly mixing with the chiral quark multiplets) and the Higgs are the lightest composite remnants of 
the strongly interacting sector, the dominant effects are expected from operators involving these fields. It is therefore 
meaningful to focus primarily on the leading dimension five contributions.^ It turns out that they can always be 
parametrized in a way that preserves the form of gauge interactions of the renormalizable theory. Therefore the only 
way to approach and constrain such terms is by studying their impact on Higgs phenomenology, which is the main 
topic of the present work.^ 

The paper is structured as follows. In Sec. II we present some general considerations of models with vector-like 
quarks including generic flavor, electroweak and Higgs constraints on their interactions, stemming from renormalizable 
as well as leading dimension five non-renormalizable terms in the effective Lagrangian. Then in Sees. HI - V we consider 
three specific model examples and analyze their viability in light of the derived direct and indirect constraints. For 
simplicity we will consider examples where vector-like quarks appear in existing weak representations of SM chiral 
quarks, thus allowing for kinetic mixing with chiral quark multiplets, with interesting phenomenological consequences. 
We summarize our conclusions in Sec. VI. Several supporting derivations and analyses of experimental constraints 
have been relegated to the appendices. 

II. GENERAL CONSIDERATIONS 

A. Renormalizable models 

Since all vector-like quark models under consideration only contain colored fermions in SM gauge representations 
and charges, we can start by considering the mass matrices of the up- and down-type quarks in the weak (chiral) 
eigenbasis 

- /:ma.s = "IXJf < + rflXd'4 + h-C. , (3) 

where the indices z, j run over all dynamical quark flavors (including new vector-like generations). The mass matrices 
■Mu,d can be diagonalized via bi-unitary rotations as A^M.d,diag = U^' -Mu.dU^' . Consequently, the gauge and Higgs 
interactions of physical quarks in the mass eigenbasis can be written in the general form (c.f. [8]) 

^^ = -^i^i^'^^PLd' + V^u'rPRdnW+ + h.c. , (4) 

Cz = -^ {xr^u'j^PLU^ ~ Xf^d'^PLd^ + Y^'^^u'rPRU^ - Y.^^d^rPRd' - 2sl,J^^) Z^ , (5) 

£,f = -{Xr^ - Y;;)'^u^Pnu^h - {Xf^ - Yf^)'^d^PRd^h + h.c. , (6) 

where Pr^l = (1 i 75)/2, g — 2mw /v ~ 0.65 is the weak coupling, while sw — ^0.23 and cw = a/1 — s^ are the 
sine and cosine of the weak angle, respectively. Jgj^j = (2u*7'^m* — d'^j^d^)/3 is the EM quark current. The flavor 
matrices V^'^, X"'"* and Y'^''^ are all given in terms of U^'r, in particular we can write V^j = {Ul)*f.{U^)ik, where the 
repeated index runs over all left-handed weak doublets, and V^^ = {UR)%{UR)ik, where the repeated index runs over 
all right-handed weak doublets. Then the (hermitian) flavor matrices entering neutral current and Higgs interactions 
are given simply by X" = V^V^^ , X'^ = V^^V^, F" = V'^V^^ and X'^ = V^^V^. Thus, non-standard Higgs 
interactions in such renormalizable models with extra quarks are necessarily constrained by charged and neutral weak 
currents among the known three generations of quarks. For example, the departures of V'" from 3x3 unitary matrix 
and the appearance of a non-zero V^ are constrained by precisely measured tree level charged current processes. 
For example, Ej=d.s,J^»j I' = 1 - Af < 1 (for i = u,c,t) and J2j=u,c,t \^ji? = 1 - Af < 1 (for i = d,s,b) are 
constrained in absence of V^ as Aj; < 0.001 [9], A^ < 0.052 [10], Af < 0.13 (see Appendix A for details) , A^ < 0.01, 
Af < 0.08 [10] and Af<l- [V^f [^ < 0.15 (see Appendix A for detafls). Note that in models with no extra up-type 
(down-type) quarks, Af = dXf^ (Af = 6X:^), where SX^^''^ = 1 - X'^f. The entries of V'^ on the other hand, 
are also constrained at the tree-level by searches for right-handed charged currents (c.f. [11] for a recent analysis). 



^ For a recent analysis of dimension six operator effects in composite Higgs scenarios without dynamical vector-like fermions see [6] . 
^ For recent related studies in the context of explicit composite Higgs model realisations see [7] . 



Unfortunately, without information on the matrix elements involving also extra quarks present in the model beyond 
the known SM generations, these cannot be directly related to Z and Higgs couplings (y"^'') . 

In addition, one can obtain tree-level constraints on the off-diagonal entries of X'^''^ and y '"^ directly from their 
contributions to Z-mediated FCNCs of up- or down-type quarks. In all scenarios we consider, either nonstandard 
X"' 7^ 8ij or Y^^ 7^ are generated but not both. In this case the bounds on non-diagonal entries oiX'^''^ or F"''* read 
|X,"J, |y,"J < 2.1 X 10-4 [12, 56], |Xr„.,J, |r,^,,J < 0.14 (see Appendix A for details); Re(X,^J,Re(yji) < 1.4 x 10-^, 
\^dhV \Xdh\ < 4 X lO"'^ and \X'l^, \Yf^ < 1 x 10"^ (see Appendix B for details). Finally, electroweak measurements 
provide strong tree-level constraints also on the diagonal entries of X'^-'^ and F"''' corresponding to the five light quark 
flavours (see Appendix C for details). 

The main consequence of the above discussion is that in renormalizable models with additional vector-like quarks, 
Higgs couplings to the known three generations of quarks, except possibly the top (i.e. AT^", F^"), must remain SM-like, 
irrespective of the spectrum or interactions of additional heavy quarks. This is because they are rigidly related to the 
corresponding Z couplings, and thus subject to severe constraints from charged and neutral weak currents. In order 
to possibly obtain more interesting Higgs phenomenology, we are thus led to consider effects of higher-dimensional 
operators. 



B. Including non-renormalizable Higgs interactions 

Treating the SM as an effective field theory with particle content valid below a UV cut-off scale A, it is well known 
that the leading higher dimensional operators involving quark fields are of dimension six, a virtue of the chiral nature 
of weak interactions in the SM. Thus, effects of NP degrees of freedom appearing above A in low energy observables 
are suppressed by at least two powers of 1/A. On the other hand, in presence of dynamical vector-like quarks, the 
leading non-renormalizable operators can appear already at dimension five. In general, they are of the form H^HQQ 
and H^ HqQ, where q denotes the SM chiral quark multiplets. The main consequences of these new interactions are 
(i) direct di- Higgs coupling to physical quarks [xij in eq. (1)] with possible implications for the SM hierarchy problem; 
(ii) modifications of single Higgs - quark couplings [i/ij in eq. (1)] not related to weak neutral or charged currents. In 
the quark mass eigenbasis these additional contributions can generally be written as 



'-h — 



A 



UlUr + 



vdi N 






h.c. , 



(7) 



where A is the UV cut-off scale of the effective theory encompasing the SM together with a number of additional 
quark-like states. First note that the appearance of A"' ' couplings of the known three generations of quarks to 
the Higgs is a manifestation of mixing between chiral and vector-like quarks, which is in general unrelated and thus 
unconstrained by charged and neutral weak currents. On the other hand, naturalness of the hierarchical quark mass 
spectrum would require \Xf-X'^[ ju^/A^ < rairrij [14]. In order to keep our analysis as general as possible, we shall not 
impose such a condition on the parameter space of our models, although one should keep it in mind. The off-diagonal 
values of A"^''' are constrained by low energy flavour observables [15]. In the up-sector, ]A^^ ^„]u/A < 7 x 10"^ and 
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(c, u)h decay searches, respectively. Similarly, 
< 2 X 10"" and 



l^.l.>/A 



\XtL,\vlK 
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^w/A < 0.34 are constrained by D^ mixing and t - 

Ar°, Bd and Bg mixing measurements require \Xf^ ^Jv/A < 2 x 10 
1 X lO-'^, respectively.^ 

Potentially the most striking tree-level effects on Higgs phenomenology in non-renormalizable vector-like quark 
models are thus modifications of the flavor diagonal Higgs couplings to lighter quarks. In particular, in the SM the 
the total Higgs decay width is dominated by the h —>■ bb channel, the first hints of which have also been observed 
at the LHC [16, 17]. The modifications of yi,i, in eq. (1) can thus have important consequences for all experimen- 
tally observed Higgs signals. Similarly, while the h — > uu, dd, ss or /i —> cc decays are very suppressed in the SM 
and also cannot be reconstructed at the LHC due to the large QCD backgrounds, they can contribute to the to- 
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Higgs decay width in case of non-zero A^"' . In particular, defining A7 = X) f 



h~^ff 



/r 



1^, we obtain 



Xj^^^ |A,j'^'u/A — mi/v| +'}Zi=uc\-^u''^/^~ ''^i/'"\ — 10 "^^7 showing that sizable enhancement in these decay 
channels is possible (although the required values of A^"' /A would necessarily violate the corresponding quark mass 



All bounds from neutral meson mixing apply in absence of large cancellations with the tree-level Z-mediated X"''*, Y^''^ contributions. 
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naturalness conditions) and that a non-trivial constraint on X^^ can in principle be obtained from the total Higgs 
decay width. This rises a question of importance of uu — > h (or to a lesser extent dd — >■ h) production mechanism 
compared to the dominant gg -^ h mode at the LHC and Tevatron. In the zero-width approximation and at leading 
order in QCD, the ratio of the relevant hadronic cross sections in the two cases can be written as 



o'iPiP2 -» h)uu _ T{h^ uu) ^p'^p^Jt) 
cr{piP2 -^ h)gg T{h^ gg) CIIp^{t) 



(8) 



where r = m\/ s with s being the invariant collider energy squared. The relevant luminosity functions at hadronic 
{P1P2) colliders are given by 
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F,l {y)IT. i^/y) + C (y)fl' i^M 



(9) 



where f^' are the corresponding parton distribution functions (pdfs). Using the LO MSTW2008 [18] set, with the 
factorization and renormalization scales fixed to m/j ~ 125 GeV, the ratio £""/£^^ for ^/s = 7TeV and y/s = 14TeV 
is 3.9% and 2.3% respectively.** Thus, even assuming comparable h ^> uu and h ^> gg decay rates, the up-quark 
contribution to Higgs production at the LHC is below the ^ 12% theoretical uncertainties [19] of the dominant gluon 
fusion production cross section. Conversely, at the Tevatron, we find the relevant ratio Cpp/Cpp for -y/s = 1.96 TeV to 
be sizable 26%. However, given the low statistics in the gluon fusion Higgs production channel at the Tevatron [20], 
this again gives no relevant constraint on T{h — )■ uu). In the future, enhanced di-Higgs production at the LHC could 
possibly offer a competitive constraint on X^'^/A . In particular, the relevant LO hadronic cross section is given by 



(^{piP2 -^ hh) 
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and Ph = a/1 — 4m^/s. Using the same pdf parameters as above we obtain a{pp — > hh)-^ /[(|X,';'4l/0-03)(lTeV/A)]^ ~ 
5(11) fb at the 8 TeV (14 TeV) LHC, compared to the SM LO predictions [21] of aipp-^ hh)^"' = 4(16) fb, 
respectively. 

C. Impact of existing Higgs data 

Most interesting effects involving light vector-like quarks in Higgs phenomenology appear at the one-loop level. In 
general, Higgs-fermion interactions of the form (1) will contribute to gluon fusion production and di-photon decay of 
the Higgs at one loop 
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|Fi(TM/) + |i^i/2(rt)| 
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where d{ri) and C{ri) are the dimension and index of the color representation of fi, respectively, and Qi is its electric 
charge. The relevant loop functions Fiir) and Fi/2(t) can be found e.g. in ref. [22], and r^ = mf^/{Am'^). In the limit 
of large fermion mass, ^1/2(1') -> -^'1/2(0) = 4/3. In the SM, gluon fusion production, gg — > h, is dominated by the top 
quark loop with (l/2)Fi/2('''t) = 0.688, with minor contribution from the bottom quark (l/2)Fi/2('''b) = — 0.04 + i0.06. 



Using instead the NNLO MSTW2008 pdfs with the same scale choices yields ratios of 4.0% and 2.5% respectively. 



Conversely, Higgs decay to two photons in the SM is dominated by the W boson loop yielding Fi{tw) — —8.34, and 
interfering destructively with the top quark contribution of (4/3)i^i/2('''t) = 1-84. It turns out that lighter quark 
contributions to loop induced Higgs processes are negligible even if their couplings to the Higgs saturate the limits 
from A7 as discussed below. 

In order to evaluate the current constraints on modified Higgs interactions, we analyze the latest available Higgs 
data, presented in Table I. Measurements are given in terms of Higgs signal strengths normalized to SM predictions 

H-A^h — SAI nSM ' y-^' 

where A — >■ /i and h —>■ B stands for different production mode and decay channel, respectively. Experimental best-fit 
values and variances are denoted by jii and af , respectively. Experimental collaborations generally provide plots with 
separate contribution from VBF plus VH, and from ggF plus ttH production channels for a given decay channel. In 
this case, we take into account their correlation, obtaining a correlation parameter p by reproducing the plots. The 
contribution from these decay channels to the total x^ function is 

~ yP'^GF^VF (f^GFJ / VMyF-MvF/ 

where GF stands for ggF+ttH, and VF stands for VBF+VH, and index i runs over decay channels.^ If the separation 
into production modes is not provided, we use the data from different search categories for a particular decay channel, 
which generally target certain production mechanism, but does not imply 100% purity. Inclusive categories are 
dominated by ggF (~ 90%), while VBF-tagged categories can have 20% to 50% contamination from ggF. VH- and 
ttH-tagged categories are assumed to be pure. In this case, we write 

CTA^h _ ^ CTggF , ^ CTyBF , ^ CTVH , ^ CTttH ,,„^ 

-ZSAT ~^9aF-sM+^VBF-s]^+t.VH-sM +^ttH-sM^ 1^°) 

"A^h "ggF "VBF "VH "ttH 

where ^i represent contributions of the specified production mechanisms for the given category. We do not assume 
correlations here, and add each search category to x^ separately. 



x^=Erv^ ' (17) 



and the total x^ function is given by x^ = Xi + xi • 

We take into account the recent evaluation [19] of ggF production in the SM at approximate N'^LO in perturbative 
expansion, which exhibits a 17% shift with respect to the values adopted by experimental collaborations, by rescaling 
central values for signal strengths which depend on ggF production by a factor 1/(1-1-0. 17^ggF)- On the other hand, we 
have checked that the resulting slightly reduced theory error (from 14% to 12%) has negligible effect on the reported 
variances. 

In the following, we present results of our analysis. The SM is in overall very good agreement with the data. The 
associated x^ for 19 observables presented in Table I is x^g^ = 16.2 corresponding to a p-value of 0.64. Within 
the vector-like quark scenarios, all the modifications to Higgs signal strengths can be expressed in terms of four 
parameters, Rgg, R^y, Rbb and A7, where 



^h^bb I \ ybb\v 
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^^^-f^=(^E^) • (18) 

^ h^bb 



In particular. 



h^ll __ ^gg p h^ZZ,WW,TT _ Rgg h^-f^i _ R^rj h^ZZ,WW.TT _ }_ h-^bb _ ^bb ,,„x 

^^GF — p -"-77 ' ^^GF — p ' ^^VF — p MvF — p ' ^^VH — p • U-'J 

The modification of the total Higgs decay width coming from Rgg, Rbb and A7 is taken into account by writing 

f = -^ = 0.569Rbb + 0.317 + 0.085i?gg + A7 , (20) 

^ tot 



ttH contribution in GF is less then 1%. 



Decay channel 


Production mode 


Signal strength 


Comment 


ATLAS 


h^ ZZ* 


Inclusive (87%ggF) 


1.5 + 0.4 


[16, 23] 


h^bb 


VH 


-0.4+1.0 


[16] 


h^WW 


ggF+ttH 


0.79 + 0.35 


Correlation p = -0.3, [16, 24] 


VBF+VH 


1.6 + 0.8 


h^-j-j 


ggF+ttH 


1.60 + 0.44 


Correlation p = -0.4, [16, 25] 


VBF+VH 


1.80 + 0.87 


h ^>- TT 


ggF+ttH 


2.2+1.6 


Correlation p — —0.5, [16] 


VBF+VH 


-0.3+1.1 


CMS 


h^bb 


VH 


1.3 + 0.7 


[17] 


h -> WW 


0/ljet (97%ggF) 


0.76 + 0.21 


[26] 


VBF-tag(20%ggF) 


0.0 + 0.7 


[1-] 


VH 


-0.3 + 2.1 


[1-] 


h^ ZZ' 


ggF+ttH 


0.90 + 0.45 


Correlation p = -0.7, [27] 


VBF+VH 


1.0 + 2.3 


h — >■ 77 


ggF+ttH 


0.52 + 0.40 


Correlation, p = -0.5, [28] 


VBF+VH 


1.5 + 0.9 


/i — > rr 


0/ljet (80% ggF) 


0.73 + 0.51 


[29] 


VBF-tag (20% ggF) 


1.37+0.63 


[29] 


VH 


0.75 + 1.5 


[29] 



Table I; LHC Higgs data used in the analysis. 



where A7 is constrained to A7 > 0. 

We consider four different scenarios, with different choices of the fitting parameters. In each case, the best-fit point 
is determined by minimizing the x^ function. Results are presented in the {Rgg, FLyj) plane, after marginalizing over 
the other parameters. We define 68.2% (la) best-fit region to satisfy x^^^ < X^ < Xmin + 2-3, and 95.5% {2a) best-fit 



region to satisfy x 
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rnin 



2.3<x'<X: 
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6.2. 



First, we take Rgg and R^j as fitting parameters, while fixing Rtt and A7 to their SM values. This applies to 
scenarios, where the couplings of SM quarks to the Higgs (the Yukawas) are not modified, while Rgg and R^j receive 
new loop contributions. In models with vector-like quarks this corresponds to the limit of zero mixing between 
the chiral and vector-like quarks. The minimum of the x^ corresponds to a point {Rgg, R^^) = (0.71,1.30) with 
Xmin — 12.7 and p-value 0.76. Second, we take Rgg, R^^y and A7 as fitting parameters, while fixing i?{,f, = 1. This 
corresponds to scenarios where the vector-like quarks do not mix with the b quark, but possibly with the lighter 
quark generations (see [30] for a recent model example). In this case, the point {Rgg,Rjj, A-f) = (0.79,1.33,0.12) 
corresponds to the minimum of the x^, with Xmin — 12.6 and p-value 0.7. The results for the first two scenarios are 
presented on the left plot in Fig. 1. For the first scenario, la and 2a contours are represented by (darker) red and 
(lighter) orange curves, respectively. 

In the third case, we take Rgg, R.y-y and Ri,i, as fitting parameters, while fixing A7 to its SM value. This case 
applies to the most studied scenarios in the literature, where the vector-like quarks only mix with the third generation 
(c.f. [31] for a recent study). The minimum of the x^ corresponds to a point {Rgg, R-yy,Rbb) — (0.70, 1.29, 0.97) with 
Xmin = 12.7 and p-value 0.69. In the last case, we take all four parameters to fit the data, corresponding to the most 
general case of vector-like quarks mixing with all three SM generations. Here, the minimum of the x^ corresponds to 
a point {Rgg,Ryy,Rbb,A-/) = (0.73,1.31,0.74,0.2) with xLn = 12-4 and p-value 0.65. The results for the last two 
scenarios are presented on the same plot. Fig. 1 right. For the third scenario, 1 a and 2 a contours are represented by 
(darker) red and (lighter) orange curves, respectively. 

The main observation at this point is that allowing for a modification of Rbb and/or A7 (non-zero mixing of vector- 
like quarks with some of the SM chiral quarks) significantly increases the allowed range of Rgg, while it has much 
less of an effect on R^^ . This has important implications for constraining vector-like quark models using Higgs data, 
since these generically predict much larger effects in -Rgg. In the following sections we apply these general results to 
a few simplest SM extensions with a single light vector-like quark state below the effective theory cut-off A. 
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Figure 1: Left: Fit of Higgs data taking Rg 
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Rjj = R{h — >■ 77) and A7 as fitting parameters. The best fit 



point (cross), la (dark gray) and 2a (light gray) regions are shown in the {Rgg,R.yy) plane after marginalizing over A7. The 
SM reference scenario is marked with a diamond. Results of the fit for A7 fixed to its SM value are given by the darker red 
contour {la region) and the lighter orange contour {2a region). The resulting prediction from the non-renormalizable model 
with an additional up-like vector-like quark, the model with a down-like vector-like and the model with a doublet vector-like 
quark in the no-mixing and negfigible isospin breaking limit are given by the continuous-blue, dotted-red and dashed-green 
curves, respectively. General predictions from the non-renormalizable model with a doublet vector-like quark are given by 
the region presented with the dashed-black contour (see text for details). Right: Fit of Higgs data taking Rgg, R-y-y, Rt and 
A7 as fitting parameters. The best fit point (cross), la (dark gray) and 2a (light gray) regions are shown in the {Rgg,R~f-y) 
plane after marginalizing over Rb and A7. Results of the fit for A7 fixed to its SM value are given by the darker red contour 
{la region) and the lighter orange contour {2a region). The resulting prediction from the non-renormalizable model with a 
down-like vector-like and the model with a doublet vector-like quark in the no-mixing and negligible isospin breaking limit 
are given by dotted-red and dashed-green curves, respectively. General predictions from the non-renormalizable model with a 
doublet vector-like quark are given by the region presented with the dashed-black contour (see text for details). 



III. SINGLET UP-TYPE VECTOR-LIKE QUARK 
A. Renormalizable model 



As a first example, we consider the SM extended by a vector-like quark pair {Ul, Ur) in the I2/3 representation of 
the SM electroweak group. In the most general renormalizable model the quark Yukawa interactions and mass terms 
can be described by the following Lagrangian 



£[?) = T/ict^Hd?,, + 2/jf gl^ujj + vW^^HUr + MuUlUr + h.c. 



(21) 



where H = iT2H* , H = (G+, {v + h + iGo)/V2) is the SM Higgs doublet, qi the SM quark doublets and u}^ the SM 
up-type quark singlets. Note that additional kinetic mixing terms of the form Ulu^r can always be rotated away and 
reabsorbed into the definitions of yu^u- Furthermore, one can, without loss of generality, choose a weak interaction 
basis where j/„ is diagonal and real. After EW symmetry breaking (EWSB) the mass matrices for up- and down-type 
quarks are 



M, 



yuv/V2 yuv/V'2 
Mu 



Md = {ydv/V2) . 



(22) 



The weak gauge and Higgs interactions of 4 (u, c, t, u') physical up-like and 3 (d, s, h) down-like quarks in this (mass) 
eigenbasis are given by eqs. (4)-(6), where V^ = 0, V'" is a general 4x3 matrix and X''- — I3X3 . Note that in this 
model Xj" = 1 — A" and that tree level constraints on the entries of A" already severely constrain the admixture of 
U within the physical u and c quarks. In particular, we find for the 3x3 sub-matrix of A" describing Z and Higgs 
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Figure 2: Upper limit at 95% C.L. on t — u (left-handed) mixing angle as a function of the u quark mass in the model with 
an up-like vector-like quark. The gray region marks the ATLAS experimental search bound on the renormalizable model using 
the u — >■ th decay signature [33]. 



couplings to known up-type quark flavours 



|^"-I|3X3< 



0.001 2.1 X 10-4 0.14 

0.0026 0.14 

0.13 



(23) 



Loop-level u' effects provide better constraints only on the mixing of the vector-like singlet quark with the top quark. 
Neglecting the small mixing with the first two generations (effectively setting y^'^ = 0) the t — u' system can be 
described by three independent physical parameters: two quark masses {mt,mu') and a single (left-handed) mixing 
angle {Otu), which are defined as [32] 



tan(26'; 



tu) 



mtmu' 



V2vyljMu 



Ml 



Muyi 



[{yir + ivhrv/^ ' 



V2 



ml 



m„ 



.t \2 



Ml + ^livl) 



iyhf] 



(24) 
(25) 



In terms of these, X^^ — c?rj, X 



Ctustu and Xl 



s?rr , where Ctu = cos Otu and Stu = sin( 



'tu- 



Presently, the most sensitive observable to nonzero stu is the p parameter, which receives a new contribution of the 
form [32] 



Ap = 



aNc ml 2 
16775?,/ m' 



^tu 



W '"'W 



-(1 



^tu) + ^tu'' 



2c: 



tu'„ 



log(r) 



(26) 



where r = m1^,/ml and we have neglected terms of higher order in m^ f^/m1 ^, . A comparison with the experimental 

bound of Ap°^P = 4^4 x lO^'' [10] yields a constraint on stu as a function of the u' mass as shown in Fig. 2. 

While the modified top quark coupling to the Higgs boson and the presence of an additional heavy quark can in 
principle impact also loop induced Higgs decays, namely h — >■ gg, /i — > 77 and h — >■ Z7, taking into account the 
above constraints on X" these effects turn out severely suppressed making it impossible in practice to distinguish the 
renormalizable model with a singlet vector-like up type quark from the SM in single Higgs production processes. 



B. Non-renormalizable models 



Extending the above renormalizible model with the leading dimension five operators containing the light SM fields 
and Ul.r as the only dynamical degrees of freedom below a UV cut-off scale A, Yukawa interactions and mass terms 
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in eq. (21) receive corrections which result in modified interactions between up-type quarks and the Higgs. One 
can manifestly preserve the exact mass diagonalization procedure of the renormalizable model by parametrizing the 
leading non-renormalizable contributions in terms of the replacement 

M^^Af^ + e /'/^^'^l' , (27) 

plus an additional Higgs-dependent 'kinetic mixing' operator 

- Ljj - Ci UlUj^. (28) 

After EWSB, the flavor structure of gauge interactions (and the associated bounds on X") in the renormalizable model 
is preserved and only the Higgs interactions in the mass eigenbasis receive new contributions of the form (7) where 
X'^' = (leading to Rtt = 1) while X'" = [/|;.[(0,0), (ci,C2)].C/]^^ Interestingly, even though X"' has no observable 
effects on charged current interactions of quarks, one can derive an indirect bound on the diagonal entries of X"' from 

CKM unitarity. Following its definition, X,"/ = {U^)z4 (ci{UR}*k + C2{U^)*4\ we note that |X,"/|2 is proportional 

to |(C/]^)j4p = 1 — Xj", multiplied by at most 0(1) coefficients in the effective field theory expansion c\,C2 and the 
unitary rotation If^. Furthermore, the general identity |(C/]j)i4p = {mi/Mij)'^\{U^)i4\'^ implies that C2 contributions 
to u, c interactions are severely suppressed. Consequently, CKM unitarity constraints on X^„ ^^ yield indirect bounds 
on the diagonal elements \X^'^\ < 0.03max(c\) and \X^^\ < 0.2max(c5^). Note that sizable contributions to A7 are 
well consistent with these indirect CKM unitarity constraints. 

The relevant modifications to Rgg and R-jj in the up-type singlet scenario come from the modified top quark 
coupling and the presence of the additional heavy quark in the loop. In the limit rrit^u' 3> ruh, which is a good 
approximation here, the relevant numerical expressions are given by 



|0.68rj, - 0.0401 V 0.0572 _ |-8.3 + l.Srj^'^ 



where to order 1/A and including possible small t ~ u' mixing 



2 2 

Ty ^yu^^+ y„'„' m!; = ' + '*^ (^* '*^ - ''''^^^ juH ~ '"' ^'* '*^ + '^'*^^ Am;; • ^^"^^ 

Above we have used the short- hand notation for c^jj = cosO'^jj^ s'^^j — sin^^j^, where the (right-handed) mixing angle 
9[jj is defined via tan^j^ — {mt/mu')tBn6tu- Thus, h ^ "f"/ and h — )■ gg are highly correlated in this set-up. Note 
that at the one loop level and in the large rrit^u' limit, contributions of renormalizable interactions of t and u' cancel 
exactly^, therefore, leading effects appear at order 0{v/A). 

The resulting predictions in the up-type singlet scenario are presented by the continuous (blue) curve in the 
{Rgg,Rjj) plain in left plot of Fig. 1. For concreteness we take |(c*,C2)|/A < 1 TeV~^ and rriu' > 640 GeV as 
suggested by direct searches [33] (see also the related discussion at the end of this section). Also, we take the t — u' 
mixing angle to be within the 95% C.L. experimental bound discussed above {stu ^ 0.2). From the fit to Higgs data 
we obtain the preferred parameter ranges for ry at 68% C.L. of Vy — 0.86_QQg (when marginalizing over A7) and 
Vy = 0.87 ± 0.08 (when fixing A7 — 0). Note that in both cases the fit slightly prefers ry < 1. 

Finally, turning to the naturalness condition in (2), for the case of a single vector-like up-type singlet quark mixing 
with the top it reads 

m^,cl^ + ml,s',^ = j[m,s,ui~c{c',^ + C2s[^)+m^,c,u{c[s[u+C2c',^)]+Oil/A'). (31) 

In the zero-mixing limit this leads to the prediction r^ = 1 — {mt/mu')'^, independent of the cut-off scale A. Interest- 
ingly, present Higgs data (exhibiting a preference for r^ < 1) are perfectly consistent with the naturalness condition. 
On the other hand, the Higgs fit results can also be interpreted in this context as imposing an indirect bound on the 
u' mass of to„/ > 360 GeV at 95% C.L.'^ 



® Deviations from the large mass limit, as well as higher order perturbative corrections can upset this cancellation. However a recent 
study of gluon fusion production in the renormalizable model with a singlet vector-like top partner at NNLO in QCD [32] has found 
such effects to be tiny, only a few percent for maximal mixing. 

'' For a comparison with the situation after the first Higgs data see [34]. 
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It is instructive to compare the above constraint to results of direct experimental searches for up-type singlet vector- 
like quarks. Interestingly, the most severe bound on u' in the renormalizable model and assuming dominant but small u' 
mixing with the top, m^' > 640 GeV [33] is given by the ATLAS experimental search using the u' — >■ th decay signature. 
In the non-renormalizable model the relevant couplings are given by ym' = StuCtu'^^u' /v + Stu{s^^c\ -\- c^^C2)v / K and 
Vu't = StuCtu'mt / V — ctu{.c^uc\ — s^jjC2)v / K. It is then easy to check that compared to u' — >■ tZ and u' — ^ bW rates, the 
1/A corrections can in principle enhance the B{u' — >■ th) in the small stu limit (in the extreme case B{u' — > th) — 1 
the present bound is then strengthened to m„' > 850 GeV [33]) but cannot reduce it significantly below its value in 
the renormalizable model. However, if u' does not dominantly decay to third generation quarks (but instead to first 
two quark generations), the current direct search constraints are relaxed dramatically (c.f. [36]) and rriu' — 300 GeV 
becomes a possibility. In summary, the Higgs fit already provides an interesting complementary constraint on scenarios 
with an up-type singlet vector-like quark cancelling the top-loop quadratic divergence to the Higgs mass. Although 
it is at present only marginally competitive with existing direct search bounds, it is far less sensitive to the hierarchy 
of mixings with the known three generations of up-type quarks (provided they are small) . 

IV. SINGLET DOWN-LIKE VECTOR-LIKE QUARK 

A. Renormalizable model 

Next we consider a SM extension with a vector-like quark pair (D]^,Dj^) in the l_i/3 electroweak representation. 
The most general renormalizable Lagrangian now contains the Yukawa and mass terms 

- rg) = i^'qlHd^ + y'^qiHu^ + vWlHDu + MdDlDr + h.c. . (32) 

The mass matrices of up- and down-type quarks after EWSB have the form (22) with the replacement u ^ d and 
U ^ D. In the mass-eigenbasis of 4 {d,s^h,d') physical down-type and 3 {u,c,t) up-type quarks the weak gauge 
and Higgs interactions are controlled by the general 3x4 matrix V^^ (again V^ = 0) defined as before leading to 
A'" = Isxs- On the other hand, now the entries of the hermitian matrix A'^ are experimentally severily constrained 
by their tree-level contributions to CKM non-unitarity and FCNCs in the down-quark sector and already preclude 
significant mixing of the vector-like down-type singlet quark with any of the SM quark generations 



|A'^-I|3x3< 



0.004 1.4 X 10-5 4 X 10-4 
0.006 0.001 
0.0057 



(33) 



We immediately observe that Higgs phenomenology in the renormalizable down-type singlet model is again indistin- 
guishable from the SM. In particular, considering only the dominant effects due to b—d' mixing and thus parametrizing 
-^bb — ^bD^ -^bd' = CbDSbD and AJ,^, — sl^, where c^^, -I- s^^ — 1, experimental constraints indicate Sbo — 0.05(4) 
(see Appendix C). This leads to maximum allowed relative deviations from SM predictions for the decay channels 
h ^ bb, h ^ gg and /i ^ 77 of 0.4%, 0.5% and -0.02%, respectively. 

B. Non-renormalizable models 

The leading higher dimensional modifications of Higgs interactions can again be most conveniently parametrized 
via the replacement 

Mz,^Mz, + c /'/^^'^'' , (34) 

plus an additional Higgs-dependent 'kinetic mixing' operator 

yielding new Higgs interactions in the mass eigenbasis of the form (7), where now A"' — and X'^' = 
C/f [(0,0),(cl,C2)].t/^^ Constraints on \{Ul)4t\'^ = 1-Xf^ lead to the following indirect bounds |A;^^| < 0.06max(cl), 
\^fs\ ~ 0.08max(c^) (both allowing for sizeable modifications of A7) and [A^^'l < 0.13max(c^). In fact, while the 
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Z ^ bb anomaly cannot be fully resolved in this model, the data prefers non-zero b — d! mixing with SbD = 0.05(4) 
(assuming negligible d! mixing with first two generations). This is enough to allow for 0(1) modification of ybb aud 
thus i?55 at order 1/A. In particular neglecting also the mb/Mjj suppressed right-handed h — d' mixing one can write 



Vbb 



rUb 



SbDCi 



'"b /, , b ^^ \ 



A 



1 -I- SbDCi 



rribh.) 



(36) 



On the other hand, b — d' mixing has negligible effects on the modifications to giuon fusion production and Higgs 
decay to two photons 



R. 



gg 



0.0572 + 10.65 -F 0.67yd/, 
(L652 



R 



~6.5 + 0.45y, 



d'd' 



■77 



-6.51 



(37) 



where C2 contributions to yd'd' dominate as yd'd' ~ —CbBCiv^jlymd' ■ The resulting predictions from the non- 
renormalizable model with a singlet down-type vector-like quark are presented by the red-dotted curves in the 
{Rgg,R-y-y) plaue lu Fig. 1. Again we have used |c2|/A < ITeV^^ and m^' > 350 GeV as suggested by direct 
searches (see the related discussion below). Allowing for Rbb ^ 1 (right plot), the preferred parameter regions for 
yd'd' at 68% C.L. are yd'd' = — 0.16lo;j4 (when marginalizing over A7) and yd'd' = —0.17^0^^3 (when fixing A7 = 0). 
Consequently, current Higgs data are not yet very constraining in this context. On the other hand, in absence of 
significant d' mixing with lighter quarks (for Rbb — 1 and A7 = in left plot of Fig. 1), the Higgs data already give 
an interesting constraint (as discussed below) on yd'd' — —0.12 ± 0.08. 

Finally, turning to the naturalness condition in (2), for the case of a single vector-like down- type singlet quark it 
reads in the small h — d' mixing limit 



-^ - CbDC2-j^ + 0(1/ A ) , 



(38) 



or equivalently yd'd' 



''bD 



{mt/md')'^ again independent of the cut-off scale A. Present Higgs data then provide 



an indirect constraint on the d' mass, which reads nid' > 330 GeV in the zero b ~ d' mixing case and grows stronger 
for non-zero Sbo- This is to be compared to direct experimental searches [35], which yield m^' > 480 GeV for the 
renormalizable down-type singlet model dominantly mixing with the b . In this case however, the direct constraint is 
dominated by the d' — >■ Wt decay signature. Enhancing the d' — > bh rate in the small b — d' mixing limit through 
the coupling yd'b — —CbDCb£)c\v/A can thus naturally relax it to irid' ^ 350 GeV; dominant (but small) mixing with 
the first two generations possibly even further. In light of this, the Higgs data already provide a complementary and 
competitive handle on such models. 



V. DOUBLET VECTOR-LIKE QUARK 
A. Renormalizable model 



As a final example, we consider the SM extended by a vector- like pair {Ql, Qr) in the 2i/g electroweak represen- 
tation. The most general renormalizable Lagrangian now contains the Yukawa and mass terms 



,(0) 



- C'^' = y'd'qlHd'j, + f^qlHu'j, + yhQLHd'j, + yhQLHu^R + MqQlQr + h.c. 
The mass matrices of both up- and down-like quarks after EWSB now have the form 



(39) 



Mu = 



yuv/V2 

yuv/V2 Mq 



Md = 



Vdv/V^ 
yDv/V2 Mq 



(40) 



In the quark mass eigenbasis the weak gauge and Higgs interactions of 4 (u, c, t, u') physical up-like and 4 (d, s, b, d') 
down- like quarks are governed by two 4x4 matrices, a unitary V^ and a non- unitary V^. Consequently AT"' — 6ij, 
while K"' are hermitian and constrained as 



|>^"|3X3< 



0.11 2.1 X 10-4 0.14 
0.018 0.14 



|>^'|3X3 < 



0.1 1.4 X 10-5 4 X 10-" 
0.21 0.001 

0.03 



(41) 
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Figure 3: Allowed region at 95% C.L. for u' — d! mass splitting as a function of the u quark mass in models with a doublet 
vector-like quark. The narrowest purple bands apply to the renormalizable model, the middle orange band stands for the 
non-renormalizable model in the zero t — u' mixing limit, while the broadest green band is for the non-renormalizable model 
with non-zero t ~ u' mixing effects allowed by Z — >■ 6b data. The gray area marks the ATLAS experimental search bound on 
the renormalizable model using the u — > th decay signature [-33]. 



Due to such severe experimental bounds on the mixing of vector-like doublet components with the first two quark 
generations, and also with the b quark, the dominant effect on Higgs phenomenology could possibly come from the 
mixing in the top sector (via induced Y^"), which remains unconstrained at the tree-level. However, as shown in [32] 
the (right-handed) mixing angles in the top {t — u') and bottom {h — d!) quark sectors are related via the mass splitting 
between the two extra quark states u' and d! as 



4[1 - 5^(1 - rla')] = ml,[l - ^^^(l - ' 



tu' J 



where Vij = rai/'Jiij . Left-handed and right-handed mixing angles are now related through tan^^ 



(42) 

'ij' tanOij. At 

the one-loop level, the u' — d' mass splitting (Amg = m„' — m^') is constrained from EW precision measurements. 
In particular, Z ^ bb observables constrain the b — d' and t — u' mixing angles as shown in Fig. 5 (see appendix C 
for details). Together with a constraint from the p parameter, this gives the bound on Attiq as shown in Fig. 3 
(the narrowest purple bands).® Taking all this into account (in particular discarding the fine-tuned solution for large 
negative Attiq), we find (in accordance with [32]) that the vector-like quark doublet with renormalizable couplings 
has unobservable effects in single Higgs production and decay processes. 



B. Non-renormalizable models 

At the non-renormalizable level, the doublet vector-like quark model offers a somewhat richer structure than 
the singlet examples. Namely, the introduction of dimension five operators allows to shift the vector-like mass 
independently for both isospin components of Q via the insertion of an iso-triplet combination of Higgs fields 



MQQRQL-^MQQRQL + '^{vy2-\H\^)QRQL ' ^' 



Qr{HH^ - HH^)Ql ■ 



A' ' ' ' '^'^^^ A 
Similarly, one can now introduce two new operators (via iso-singlet and iso-triplet Higgs field insertions) 



£(1) = ('^l )\..2 /o ,z.,2^^_„, , (cD' 



A 



-ivy2 -\H\')QRql 



A 



Qr{HH^ - HH^)qi . 



(43) 



(44) 



* Note that we do not find a bound as strong 
quaUtativcly unchanged. 



as reported in [32], The resulting implications for Higgs phenomenology remain however 
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The two isospin breaking corrections (proportional to c]~2) now necessarily induce corrections to quark masses and 
mixings. In particular, the resulting changes to Mu,d can be parametrized as 

iyu.,y^ - i,:j- - (...)- §^ T ^^%#^ , (45) 

{yu,Dy -^ iyu,Dy = [yu.oyjY^ ± ' . ' , (46) 

Mu,D ZAAlu^D 

Mq -> Mu,D , (47) 

where Mu,d = J M^ j^ + u4((c-)i(cj;)*)/4A2 and Mu^d = Mq ± v^c^^.!^^- The only observable effect of these shifts 
is the breaking of correlation between the masses and mixing angles in the up- and down-type quark sectors , such 
that Atoq = rfiu' — rud' becomes an independent free parameter, given in the zero-mixing limit (when y'jj u ^ 0) 
solely by Amq — w^Cj^/A. At the one-loop level it will affect the p parameter as 

where we have only kept the leading Atuq dependence. The resulting constraint is shown in Fig. 3 (in middle orange 
band). However, if we also include non-zero t — u' mixing effects, marginalizing over the allowed values of Stu from 
Z ^ bb data we obtain a much weaker bound on Amg shown in the uppermost (green) band in Fig. 3. In our 
numerical evaluation we employ the full one-loop formula for Ap, which can be found in Appendix D. We thus 
conclude that after including the contributions of leading higher dimensional operators, the isospin components of a 
TeV scale vector-like quark doublet can be split by as much as 30%. Although this has no observable consequences 
for Higgs phenomenology, it can have profound implications for direct m' searches if the u' — > d'W decay channel 
becomes kinematically allowed. 

The new Higgs interactions in the mass eigenbasis are again of the form (7), where now AT"'''' = U]^' .[(0, 0), (c^ ± 
cr,c+ ± c^)]*.C/^''^ or explicitly {X^'^, = {U^%[{ct ± c{?{Ul% + (c+ ± c^){Ul%]. They are thus 
constrained indirectly by bounds on Y"'^, since now \{U'^ )n\^ — Y^' . Taking into account also the relation 
|(C^l' )j4p = i'mi/Mu^D)'^\{U^' )i4pi we can safely neglect cf contributions and obtain X^'^ < 0.35max[(c^ -I- cj~)*], 
X^"J < Q.13m&x[{ct +c^y], Xf^ < 0.30max[(c+ - c^)*], Xf, < 0.40max[(cf - q)*] and X^l < 0.11 max[(c:;^ - q)*]. 
Significant effects in A7 and Rtb are thus possible and can be used to constrain the diagonal entries of X"'"^' . 

Additional interesting effects again appear in loop induced Higgs processes. Higgs decays to pairs of gluons or 
photons are modified by additional heavy particles in the loop 

_ |0.68(r,+r^)- 0.0401' + 0.057' _ |-8.3 + 1.8r, + 0.45r^|' 

where to order 1/A and including t — u' mixing 

V V 

Tx = ytt i-yu'u' — 

2 2 

= l^stu {{ct + c^yc',^ - (4 + <h)-^tu) -^^ - <^tu {{ct + c^Ys'.u + (4 + c-)c;^) £-- , (50) 

and 

ry = yd'd' — = -T {4-c:^). (51) 

md' Anid' 

Taking into account the bound on t — u' mixing, a good approximation for r^ = 1 + Stu ^m Km ■ There 

is no correlation between Rgg and R^^ in general, unless one imposes additional constraints on the parameters. 

Therefore, the resulting predictions from the non-renormalizable model with a doublet vector-like quark are given 
by a region (dashcd-black contour) in the {Rgg,Rj^) plane of Fig. 1 left, in the case when modification of A7 is 
allowed and Rbb — 1, and on Fig. 1 right, when sizable modification of R^b is allowed as well.^ We have assumed 



^ As in the case of the non-renormaUzable model with a down-like vector-like quark, ytb can receive 0(1) modifications even for small 
b — d! mixing. 
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ki 2~| 1^ — 1 TeV~ and rUqi > 790 GeV as suggested by direct searches [33]. Also, we take the t — u' mixing angle 
to be within the 95% C.L. experimental bound discussed above {stu ^ 0.35). 

For small mixing and negligible isospin breaking, masses of u' and d' quarks are degenerate and Vy — r^ — I = 
~ Am~''2 • Allowing for modification of A7 and fixing Rbi, — 1, the resulting predictions in this scenario are presented 
by the green-dashed curve in the {Rgg, R^^) plane in Fig. 1 left. Allowing for sizable modification of Rbb as well, 
the predictions are presented by the green-dashed curve in the (Rgg, R-y-y) plane in Fig. 1 right. Again we have 
assumed |cj| /A < ITeV^^ and mq> > 790 GeV. In particular, the preferred parameter regions for Vy at 68% C.L. 
arc Ty = —0.09^006 (marginahzing over both Rbb and A7), ry — — 0.09^Q'og (marginalising over Ri,b but fixing A7), 
Ty = —0.07^004 (fixing Ri,b = 1 and marginalizing over A7) and finally ry — —0.06 ± 0.04 (fixing both Rbt and A7 
to their SM values) . In this scenario, the Higgs mass naturalness condition reads 

4 =24"^+0(l/A^), (52) 

or equivalently Vy — —{mt/V2mq')'^. This condition allows to put an indirect bound on the mass of q' to be 390 GeV 
at 95% C.L. in the case of fixed Rbb- 

Turning to direct searches, the most severe bound on u' in the renormalizable model with a doublet vector-like 
quark assuming dominant but small mixing with the third generation, m„/ > 790 GeV [33] is given by the ATLAS 
experimental search using the u' — > th decay signature. In the non-renormalizable model the relevant couplings 
are given by ytu' = StuCtum^'/v + Stu{s'tij{cf + cj")* + c^jj{c^ + c^))«/A and y^'t = StuCtumt/v - Ctu{c'tu{c^ + 
c^)* — s^u{c2 + c^))z;/A. It is then easy to check that again compared to u' — > tZ and u' — > bW rates, the 1/A 
corrections can in principle enhance B{u' — >■ th) in the small Stu limit but cannot reduce it significantly below its 
value in the renormalizable model. However, the significant u' — d! mass splitting allowed by present data when 
including dimension five contributions, reopens the possibility that the dominant decay channel of u' is actually 
u' — > d'W , in which case the existing direct search constraints are considerably relaxed and dominated by searches 
for the lighter isospin component via d' — >■ tW and d' — > bh decay signatures (in the case of dominant mixing with the 
third generation). Consequently, in such scenarios u'{d') masses as low as rn^i(^^i^ > 400 GeV (300 GeV) could still be 
viable. 

VI. CONCLUSIONS 

We have systematically investigated the impact of dynamical vector-like quarks, accommodated within SM gauge 
representations and charges, on Higgs physics. In particular, we have considered the weak singlet up-type, singlet 
down-type and doublet vector-like quarks, potentially mixing with all three known generations of chiral quarks, and 
have updated the most relevant constraints on such scenarios from low energy flavour phenomenology and electro- weak 
precision measurements. Within the renormalizable SM extended by additional vector-like quarks, we have shown 
generally that Higgs couplings to the known three generations of quarks need to remain SM-like regardless of the 
extra quark masses. This feature is a consequence of the fact that precision flavour and electro- weak observables are 
affected by the mixing of vector-like and chiral quarks, some of them already at the tree level. Consequently, Higgs 
decay widths to light quark pairs, gg and 77 cannot deviate significantly from their SM predictions. 

A singular feature of models with vector-like fermions is that non-renormalizable contributions sensitive to physics 
at the cut-off scale of the effective low energy theory appear already at dimension five. The inclusion of such higher- 
dimensional operators is essential in models that aim to cancel dominant quadratic divergences to the Higgs boson 
mass coming from top quark loops with new fermionic contributions. Contrary to the renormalizable models, they 
also predict interesting effects in Higgs phenomenology. We have investigated such contributions for all three types of 
additional vector-like quarks (see Fig. 1). The most important consequences for Higgs physics are possible significant 
enhancements in Higgs decay rates to pairs of light {uu, dd, ss, cc) quarks, which may still account for a significant 
fraction of the total Higgs decay width. This is possible due to a simultaneous strong modification of the Higgs 
production by gluon fusion gg ^ h (up to 50% deviation from SM predictions, while the /i — >■ 77 decay width can 
receive modifications up to 10%). Such a possibility could thus be tested by future more accurate determinations of 
vector boson fusion and {W, Z or ft) associated Higgs production. In addition, a possible modification of the Higgs 
coupling to b quarks in models containing down-type vector-like quarks will be probed by future precise measurements 
of the h ^ bb decay. 

Interestingly, current Higgs data are perfectly consistent with (and even exhibit a slight preference for) the possi- 
bility that vector-like quarks contribute to the cancellation of the top loop quadratic divergence to the Higgs mass. 
Conversely in some cases, a fit to existing Higgs measurements under such an assumption already offers competitive 
and robust constraints on vector-like quark masses in comparison with results of direct experimental searches which 
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need to rely on particular decay signatures. In the example of a weak doublet of vector-like quarks, we have shown 
that dimension five contributions allow to relax the stringent bounds on the mass splitting between the two isospin 
states. Consequently the decay of the heavier doublet component to the lighter one with the emission of a M^ boson 
may become kinematically allowed, affecting the relevant experimental signatures. More generally in presence of 
dimension five contributions, vector-like quark decay widths are naturally dominated by decay channels involving the 
Higgs. Future dedicated experimental searches for such particular signatures (as exemplified in [33], see also [37]) 
could thus shed light on the relevance of vector-like quarks in the solution to the SM hierarchy problem. 
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Appendix A: CKM non-unitarity and Z mediated FCNCs in top quark production and decays 

In absence of right-handed charged currents, experimental constraints on X]i=d s h \^tiV = 1 + A" < 1 and also 
[Vj^ 1^ can be obtained from the measurements of the (i-channel) single top production cross-section (at) at the LHC 
and the fraction of top decays to Wh pairs in tl production (i?). Assuming t — > Wq channels dominate the top 
decay width, to a very good approximation R is given by i? ~ 1^*6 P/(Si=rf s b l^tf P) ^^^'^ '^^ '^^^ ^^^ ^^^ recent CMS 
result R'^^'P = 0.98 ± 0.04 [38]. Note that this measurement alone requires that \V^^\ > |Vt^|, \V^^\. The relevant 
t-channel single top production cross-section can then be written as at — af^\V^j^\'^R, where the SM prediction of 
af^ = 64.6^2^0 pb [39] is obtained with \VtJ^\ = 1 and i? = 1. We compare this to the weighted average of the recent 
ATLAS [40], and CMS [41] results (u™'' = 68.5 ± 5.8 pb ). Performing a x^ fit of the two experimental quantities in 
terms of |Vj^|^ and J2i=d s b l^tf P {a.nd taking into account theoretical constraints J2i=d s b l^fP — ^ ^'^'i \^u\ > 0) 
we obtain the 95% C.L. lower bounds of 

|V',[|2>0.85, Y. |Vii^l'>0.87. (Al) 

i—d.s,b 



In models with up-type weak singlet vector-like quarks, \Xy,t\ and \Xct\ will contribute to FCNC top decays. 
Combined with at and R, the experimental bounds on B{t — )■ Zq) can then be used to constrain -^/iXuf [^ + \Xct\'^ = 
|A"t„^tc|- The presence of these new decay channels in principle also needs to be accommodated in at and R by writing 



<rsM\Vti\'{^+Pwz^-^^) , (A2) 



where 

1 (2Af|-hm?)( 



Mi 



Pwz^Ty ; ^77T^' (^^2) 



takes into account the dominant phase-space difference in i — ?► Wq and t ^>- Zq decays. In addition, searches for 

/"jf I = 1. In presence of [V^f | 



t —>■ Zq typically assume B{t -^ Wb)+B{t —>■ Zq) = 1 and | V^^l = 1. In presence of \V^^\ < 1, the experimental results 



should instead be compared to 



rL\2 N -1 



Including the recent ATLAS result B{Z -^ tq)''''^ < 0.73% [42] in our fit, wc first observe that the presence of FCNCs 
has no observable effect on the results in eq. (A2). On the other hand, we can obtain an upper bound on lATttitcj < 0.14 
at 95% C.L. . 

Although the presence of right-handed charged currents complicates the analysis of top production and decays, a 
robust bound on [ytu^tcl and also |Vj^| can nonetheless be obtained. To a good approximation namely in this case 

R = l^^'I' + l^^^l' f A5) 
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and thus experimentally l^^tP + |V(f P > IV^^P, l^^f P, iV^^lM^^f P ■ In addition, the presence of V^f wiU affect single 
top production as 

^t = ^SMd^iP + >^R\Vi\) (^1 + PWZ ^JJ^^^"^^^, ^ , (A6) 

where k^j ~ 0.92 [43] . Finally, V^f contributes to the positive W helicity fraction {J-+) in top decays as 

x-, ^ l^tbl x-SM I R \Kb\ "^^ /A7^ 

where J-^ = 0.0017(1) [44], x = [mw/'mt)^ and ?]^ — 0.93 parametrizes NLO QCD corrections [45]. Using the 
recent determination oi J^'^^ = 0.01(5) by ATLAS [40] we obtain |Vjf| < 0.2|Vj^|. Thus we may conservatively use 
the results of the previous paragraphs also to constrain |Vj^| > 0.85 and |yt„,tc| < 0.14. 

Appendix B: Bounds on down- type quark mixing witli vector-like w^eak singlets and doublets from rare K 

and Bq processes 

Model independent bounds on the off-diagonal entries of X'^,Y'^ in models with additional vector-like down- type 
weak singlet and weak doublet quarks respectively can be obtained from their tree-level (Z-mediated) contributions 
to FCNCs involving down-type quarks. For example, a bound on X^^,Yf^ can be extracted from the Kl — >■ M"*"/^" 
decay. We use a conservative estimate for the pure short distance branching fraction B{Kl — > fJ'~^fJ'~TsD < 2. 5 x 10~^, 
obtained using dispersive techniques [47], as a la upper bound. Neglecting the much smaller SM contributions, the 
X^^ contribution can be written as 



/^2 / 4j^2 

B{Kl ^ ^i+t^-)'' = Y^fKrr^KTK.mlJl- —^Kc{X,af ■ (Bl) 

Using the inputs for Gp, masses and lifetimes from [10] and also /x — 155.37(34) MeV [9] we obtain Re(X5d) < 1 x 10^^ 
(the same result apphes also to Ysd) ■ Since much stricter bounds are expected on lm.{Xsd) (J^^iXsd)) from the precise 
knowledge of ck, we interpret the above values as conservative constraints also on the moduli of Xsd and Ysd- 

In the Bd sector, Xm contributes at the tree-level both to Bd — >■ /^"""M"; as well as in B^ — B^ mixing. Neglecting 
Xjjd loop level corrections due to CKM non-unitarity and dynamical vector-like quarks running in the loop, the 
Bd -^ M^/^~ branching fraction can be written as 



q2 

B{Bd -^ ^J-~^^J.~) = -^fl^mB^TB^m 




y 2 

\t ^SM , ^bd 
^bd'^dB=l "I 1^ 



(B2) 



with A^^ = ^tdYtb* ttie relevant CKM combination and the SM Wilson coefficient given by 

CllU = -j^VvYo (4-) ■ (B3) 

Here Iq is the relevant Inami-Lim loop function [48] 

A — X 3x 



Yo{x) 



l~x (1 



In:; 



(B4) 



while rjY — 1.01 [49] parametrizes higher order QCD corrections. For the values of SM input parameters (in particular 
ctcm, sw and Wit) we follow the prescription of [49] , and we use Jb = 190.6 ± 4.7 [50] . In the B° — B^ system, the two 
most relevant observables are the mass-difference between the two Bd mass eigenstates {Anid) and the CP violating 
phase in the mixing {/3d). Since the corresponding width difference (AF^) is small [AFdl ^ [Amdl, one can write 

Am,^2|Mf,|, sin2/3, = ^^|^, where Mf, ^ ^tubJ^b^Bb, I^Xld'cj^=, + ^) ■ (B5) 



with the SM Wilson coefRcient given by 
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V271 



VbS 



'w 



(B6) 



Here again S is the relevant Inami-Lim loop function 



S(.) = I 



1 3 1 - 3a; 



3a;2 



2(1- a;)2 (1 - , 



lux 



(B7) 



Bs. =0.0411 ±0.0075 [51] 



while rjB = 0.939 [52] parametrizes higher order QCD corrections. In addition we use fg 

In absence of CKM unitarity, we need to determine not only X) 
purpose we use the radiative rate B — >■ Xd^, which is unaffected by Xm at the tree-level. Although it receives 
non-standard contributions proportional to Xm due to CKM non-unitarity and extra down-type quarks in the loops, 
these are parametrically (loop) suppressed compared to tree-level effects in B^ — > M^A*~ and B'^ — B'^ mixing. Using 
the SM recent evaluation [53] 



but also the CKM combination A^^^. For this 



B{B ^ Xrf7)|^>i.6GeV 



\t 
■^bd 



0.0084 



1.54(26) X 10" 



(B8) 



we extract \Xld\ 



and comparing it to the experimental result B{B — >• Xdj^^^^i QQ^y = 1.41(57) x 10 ^ [54] 
S.OljQ ^ 10^"^ . Plugging this into eq. (B2) we observe that compared to the experimental 95% C.L. upper limit 
of B{Bd. — > i.i^fi^)°^^ < 9.4 X lO""'^" [55] the A^^ contribution can be safely neglected and we obtain a bound on 
\Xbd\ < 4 X lO^"*. Note that the measurements of Amf^ = 0.507(4) ps^^ and sin2/3'=''P = 0.679(20) [5G] cannot be 
used to impose a stricter constraint since the phases of A^^ and X^d can always be arranged such that cancellations 
between these contributions weaken the prospective bounds above the one by B{Bd — >■ fi'^^J,^)°^^ . 

Finally in the Bs sector we can proceed similarly, with obvious replacements A^^ — > A^^ and Bd —>■ Bg- We 
again employ B — >• X^-^ to extract \X\J\ (modulo loop-suppressed X^g effects) as \\\g\ = 0.043(2) [10] and use 
Jb, — 227.6 ± 5.0 [50], f% Bb, — 0.0559(68) [51] for hadronic inputs. By comparing to the experimental values of 
Amf P = 17.719(43) ps"\'2/3f p = -0.1(6.1)° [56] and 6(5^ ^ M+M")'"'^ = i^-^-H) x 10~^ [55] we thus obtain a 
bound on \Xbs\ < 1 x lO"-^ , which is again dominated by the muonic Bg decay rate. 

The presence of vector-like weak doublet quarks induces right handed charged and neutral currents among SM 
quarks. The resulting flavor phenomenology is very rich and deserves a dedicated study. For our purpose however, it 
suffices to show that to a first approximation, one can actually neglect all terms coming from right handed charged 
current operators as well as the ones containing extra u', d' quarks in the loops contributing to quark FCNCs. This 
requires some knowledge of the mixing matrices, which can be expressed through the parameters of the Yukawa sector 
in terms of an expansion in ratios of mass parameters, which enables us to connect the left and right handed mixing 
matrices. However, the approximation only holds if such an expansion is justified, as we will check a posteriori. 
First recall that the up- and down-like quark mass matrices in presence of a single vector-like quark doublet can be 
written in the form (40). Diagonalizing the products A^A^^ and TW^A^ one obtains the left and right mixing matrices, 
respectively, which we write in the form (c.f. [57]) 



jjUjd 
^ LM 



T^u,d r^u.d 



■^L,R 



T\ 



(B9) 



L,R 



The mixing matrices between the vector-like and chiral quarks can most easily be obtained by starting from a basis of 
right handed chiral quarks, where yu,dy\ d^^/-^ ^^'^ both diagonal, which can be done as the right handed chiral quarks 
are isosinglets. The mixing yjj.DV j \f^ is bounded by the EW scale w, whereas the Dirac mass vtiq is experimentally 
required to be larger. Taking thus v/Mq as the expansion parameter, it can be shown that, to first order in this 
expansion, both 4x4 right-handed rotation matrices schematically take the form 



^ R. 




(BIO) 



so the right handed charged currents are suppressed with Vif ^ 0{v'^ /Mq) and Vj^ 
light three generations and the extra heavy quarks, respectively. 



0{v/Mq), where I, h stand for 
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Turning to the left handed sector, we can choose e.g. a basis of left handed quarks where yu is diagonal and real 
and use a unitary transformation to diagonalize yd- Then we proceed in a similar way as before, and we obtain 
that, to first nonvanishing order in v/Mq, K^' from (B9) equal those matrices (no corrections to unity or unitarity). 
Combining the up and down rotations, the 3x3 upper left submatrix of V^ takes the form 

V^'i = iKlUiKt)*k + iRlURi)h hJ,k ^1,2,3. (Bll) 

In order to obtain the matrix elements, one has to solve the equations /A^-'^M'^''^^U2' ~ U^' ^u d' where Dq = 
diag(rnq,mqi) are the diagonal mass matrices. In the process we require that the entries in U^' mixing chiral and 
vector-like components (R^' and S^' ) are smaller than the remaining ones {K^' and T^' ) in terms of v/Mq scaling 
(similar to the case of U^ ). Also, we assume that the corrections to the masses are small enough so that mq,mq> 
are of order v and Mq, respectively. Consequently the equations for i?^' read 

^'iVu^yLRl'' + Vu,dyloTL') = 2ml,^,,Rl'. (B12) 

With the mentioned assumptions, one can neglect the first term in the above equation, obtaining R^ — 
v^yu.dyij dI'^''^'u' d" Thus, R^ are one order higher than i?^' (see BIO), the off-diagonal elements in the fourth 
row and column are of higher order as well, which makes K'l a unitary 3x3 matrix to second order in v/Mq. 
Consequently all right handed contributions, as well as corrections in the left handed sector are doubly suppressed in 
the loops of Bs — > ^1^ ^l^ and we can deduce a conservative bound on the FCNC's by neglecting the afore mentioned 
corrections, as we have done it in the vector-like singlet quark case. In fact, with the approximations made, the upper 
bound on \Y},s\ is the same as the one on |Xf,s| in the down-type singlet case. 

Appendix C: Constraining Z ^^ qq 

The appearance of Z-mediated FCNC's is generically connected to modifications of diagonal Zff couplings. In the 
quark sector, such effects can be probed by direct measurements of the hadronic Z decay width, its heavy flavor decays, 
such as Z — ?> bb, but also at low energies by e.g. atomic parity violation (APV) measurements. In general we can 
parametrize the chiral Zqq couplings in terms of Q^j (M — L, R) as the coefhcients multiplying — ^^^g , (qm 7^ Qm ) Zf^ 
and SGIj = GIj - {Gm)sai where {Gm)sm are given in ref. [58]. In our scenarios SGl == 5X^^/2 or SGl = -6X^^/2, 
and 6G'n = Y^V^ or <5^^ = -r4/2 . 

Stringent constraints on Ar^„ and X^^ can be derived from APV measurements in ^'^'^Cs. The tree level modification 
of the Z boson couplings to first generation quarks will modify the nuclear weak charge as [59] 

6QwiZ, N) = 2(2Z + N){SGl + 5GI) + 2{Z + 2N){5Gi + SG^) , (CI) 

where the measured value deviates from the SM expectation by 1.5 ct [60] 

Qw - Qw' = SQw = 0.65(43). (C2) 

For the singlet up-like vector-like quark, this constrains (5A"„ = 0.0035(23), while for the singlet down-like vector-like 
quark SX^j^ = —0.0031(20). In the doublet case, one has two independent variables SC^ and SGJ^ at the tree level, 
and only a certain linear combination can be constrained, namely Y^^ — 1.12yj^^ = 0.0035(23). 

Additional constraints can be derived from direct Z ^ qq measurements presented in Table II. In particular, the 
total Z decay width is given by 

Fz = Finv + Flop + Fhad- (C3) 

It incorporates the decays to leptons, where Ficp = 251.7 MeV in the SM, decays into invisible particles (neutrinos), 
being Finv = 501.6 MeV in the SM, and the decays into hadrons. The hadronic width is the sum over the decays into 
all kinematically accessible quarks, that is all SM quarks but the top, Fhad = J2q ^g- The partial Z-decay width to 
light quarks is given by 

F, ^ F(Z ^ qq) ^ Nc^^ {Rt. \Gl + G'k\' + R\ \Gl - GU") + ^Iw.QCD, 



19 



Observable 


Measured value 


SM prediction 


Reference 


Rb 


0.21629(66) 


0.21474(3) 


[61, 62] 


At 


0.923(20) 


0.9347(1) 


[10, 61] 


aU 


0.0992(16) 


0.1032(6) 


[61, 62] 


Ohad [nb] 


41.541(37) 


41.477(9) 


[10, 61] 


rz[MeV] 


2495.2(2.3) 


2495.5(1.0) 


[10, 61] 


Re 


0.1721(30) 


0.17227(4) 


[10, 61] 


A, 


0.670(27) 


0.6680(4) 


[10, 61] 


A'pB 


0.0707(35) 


0.0739(5) 


[10, 61] 



Table II: Z pole observables used in the analysis and their SM predictions. Correlations between observables are neglected. 
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Figure 4: Left: Fit of Z-pole data taking (5G|_ ^ for q — u,d, c, s, b as fitting parameters. Best fit point (cross), 1 a (dark gray) 
and 2(7 (light gray) regions are shown in the (^G/j, cSG^) plane after marginalizing over the other parameters. The results of 
the fit for fixed 5G^ ^''^ = are given by the red contour (Ict region) and the orange contour (2ct region). Right: The fit of 
Z-pole data in the model with a singlet down-like vector-like quark. Best fit point (cross), la (dark gray) and 2a (light gray) 
regions are shown in {SX^f^SX^s) plane. 

with radiator factors Ry and R\ and non-factorizable radiative correction parameters ^%vi/qcd given in [.58]. The 
fractions of hadronic Z decays involving b quark pairs and c quark pairs are defined as 



Rb 



had 



Re 



T 



had 



(C4) 



respectively. The associated bottom and charm quark left-right asymmetries [Af, and A^), and forward-backward 
asymmetries {A'^p^ and A^g) can be written as 



Af 



2yr^4^ f^^ 



f 
Ok 



% + (l + 2z,)(|±|)' 






:A,Af 



(C5) 



1 



where f = b,c and Zf = m?Amz)/m\. The electron asymmetry parameter is fixed to its SM value, {Af.)sM = 0.1464. 
Another interesting quantity is the hadronic e+e" cross section at the Z pole (chad)- H can be written as 



Chad 



127r FeThad 



Ml 



^ z 



(C6) 



where Tg = 84.005 MeV and M^ — 91.1876 GeV. We perform a x^ fit of the data presented in Table II in terms 
of 5G\., as fitting parameters. The results are presented in the {5Q\, ^Ol) plane, after marginalizing over all other 



by 5.6, 4.9 and 1.3, respectively. In the SM reference scenario xIm — 1^.9, 



variables, in Fig. 4 left. At the best fit point Xmin — 3-''', and the most important observables in the fit are i?;, 
and Ofiad^ which contribute to x%m ~ Xr, 
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corresponding to a p-value of 0.06. The largest contributions to Xsm come from i?h, A'^pg and ahad and are 5.6, 5.2 
and 2.8 respectively. We also perform a x^ fit only with i5G^ and SG'j^ while putting other parameters to zero. In this 
case Xmin — '^■^ with p-value 0.25. The corresponding 1 a and 2ct regions are presented in the left Fig. 4 by red and 
orange curves, respectively. As expected, the data is mainly sensitive to bottom quark couplings leading to model 
independent bounds of SG^ = 0.002(2) and (5^^ = 0.015(6). ' 

Now, we turn to specific models. In the singlet up-type vector-like quark model, tree level modification of SOf, and 
6Q1 is possible. Taking into account the severe constraint on SX^^ = —0.0001(6) from CKM unitarity, 6Q2 has no 
significant influence in the fit. Therefore, we use the data to extract the best current constraint on SX^^ = —0.0020(13). 
The constraint comes essentially from three observables: Rb, <Jhad and F^, which are more constraining than direct 
Z ^ cc measurements Re and A'^pg. We calculate the individual contribution to Ax^ from each observable in the 
points which are ±1 u away from the best fit point. Contributions to ^'X^_„ from i?^, Ohadi ^ z ■, Re and A'^pg are 

— 1.1, —0.7, 2.5, 0.4 and —0.1, respectively, while contributions to Ax+g. from i?b, a^ad, F^, Re and A'^p^ are 1.4, 1.1, 

— 1.4, —0.2 and 0.1, respectively. 

In the singlet down-type vector-like quark model, tree level modifications of SQ'^ , SQ^ and 5Q\ are possible. 
APV puts a strong constraint on 5Xf^, so that the Z data can be fitted with only two parameters, SX^^ and 5X^^^. 
Results are presented in Fig. 4 right. The best fit point now corresponds to {SX^^^SX'^,,) = (0.0027,-0.0046) where 
y^nin — ^■'^ ■ The fit is statistically better than in the SM with a p-value of 0.2. The most relevant observable in the fit 
is i?(, and its contribution to x^gj^ — Xmin i^ 5.6, while the contribution of (Jhad is 1.0. Analyzing one-dimensional x^ 
functions for each observable, we find 1 tr preferred regions to be SXf^ — —0.005(2) and SXj^f^ — 0.0027(15). Finally, 
including SX^^ in the fit and taking into account the APV constraint, the preferred region for SXf^ is reduced to 
SXf^ — —0.002(3), while the preferred region for SX^/^ is unaffected. 

In the doublet vector-like quark model, tree level modification of all right handed couplings of light quarks with 
Z boson is possible. Therefore, we fit Z-pole observables together with the APV constraint on Y^^ — 1.12Fj^ with 
five parameters Y^^, Y^^, Y^^, Yf^ and Yj,^. Marginalizing over four parameters, we get the preferred range for the 
remaining one. Modification of charm and bottom quark couplings can be constrained to a percent level, Yjfi^ = 
—0.018(6) and Y^^ = —0.003(9). Negative value for Fj,^ is mainly driven by Rh and A''p^, which contribute to 
Xsm ~ Xmin by 5.5 and 3.5, respectively. In the case of the light quark couplings, observables given in the table can 
not distinguish between different light flavors, giving very poor constraints on one coupling after marginalizing over 
others. Therefore, we include new observables into the fit which are poorly measured but can distinguish between light 
quark flavours, namely asymmetries associated with the strange quark, with experimental values {A'')exp = 0.895(91) 
and {Apg)exp = 0.0976(114) [10]. The constraints we get are rather mild, since the experimental values for Ag and 
Apg have large experimental uncertainties. We obtain the following bounds on Y"^ = 0.035(40), Y^^ = 0.030(35) and 

Finally, in the model with one vector-like quark doublet mixing predominantly with the third generation, modifi- 
cation of Sg'pj is induced at tree level Sg'^ = — (l/2)s^^, while modification of Sg'l is induced at one-loop level 

b _ 9^ 2 ffiix,x') 
64^'*^ 

where x = m^/my^, x' = m^^ijm^ and r = rn^i/ml, and /i and /a are loop functions given in refs. [32, 63]. The 
above expression is given in the a;,x' 3> 1 limit and for small mixing angles. In our numerical evaluation we use the 
full one-loop expressions and do not assume small mixing angles, even though we note that the above approximations 
are fairly good. Using the available Z-pole data to fit Sbo and Stu parameters, we present the results in the {sbD, Stu) 
plane in Fig. 5 for fixed m„' = 800 GeV. We have checked, however, that the results are mostly insensitive to the 
precise value of the u' mass in the interesting region (500 GeV < to„' < 1500 GeV). Marginalizing over Sbo, we obtain 
\stu\ <0.35 at 95% C.L. . 

Appendix D: Doublet vector-like quark contributions to p parameter in presence of 1/A corrections 

In the renormalizable doublet vector-like quark model, the divergences in the loop calculation of the p parameter 
vanish only after imposing the connection between masses and mixing angles in the up- and down-quark sectors 
namely 

Md = Mu = Mq , (Dl) 

with 

M2 2 2 I 2 2 

D ~ (^bO'^^d' + ^bD''^b i 

M^ = c^^jj^l + sljjml . (D2) 



^91 = £-,s'tu '-^^^^^ + h{x.x') , (C7) 
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Figure 5; Fit of Z-pole data in the model with a doublet vector-like quark mixing predominantly with the third generation. 
Best fit point, la (dark gray) and 2ct (light gray) regions are shown in the {sbD,Stu) plane. 

Relation Dl no longer holds after the inclusion of non-renornializable operators. However, one can still relate the 
parameters in the up and down sector through the identities (see Sec. V) 



Mu^D = JM^^ + t;4((cr)»(cr):)/4A2 , Mu^D = Mq ± v'c^/2A. 



(D3) 



Expanding to 0(1/A), the divergences in p again cancel. Building upon known oblique corrections from vector-like 
quarks in the renormalizable model [G4] , the new physics contribution to the p parameter including leading dimension 
5 non-renormalizable operators can then be written as 



Ap. 



aNc 
16^s2^ I ,_, 



\ i—t,u',j^b,d' 



SbDCbD9ii^d',Xb) - s^jjCtjjgi{xu',xt) - gi{xt,Xb) + g,. 



(D4) 



where Xi = mf/rn^. The corresponding mixing matrices are defined as {V^)ij = {U^)ik{U^)*f^ and {V^')i 



iU:^UU^)h^ where 



fiR _ / Ctu,bD Stu.bD \ fjL 

'-'u.d — \ 1 '-^u.o 



-StU,bD Ctu,bD 



^tUMD ^tU,bD 
^^tU,bD ^tU,bD 



The relevant loop functions are given by 



2xia;2 , xi 
gi{xi,X2) = xi +X2 In — , 

Xi ~ X2 X2 

g2[xi,X2) = 2Va;ia;2 In 2 , 

\Xi - X2 X2 



(D5) 

(D6) 
(D7) 



while the new term which incorporates explicit effects due to 0{1/A) non-rcnormalizable mass-splitting operators 
reads 



5nr = 4:6 < sljjXb In sljjXt In ^ -f xt/ In 

[^ Xfl' Xy^' X^' 

with 5 = {Mjj — M£))/Mu and Afj/^D given in D2 (note that xu refers to M^/m'^). 



(D8) 
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